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Abstract 

New fermions, strongly coupled to the Standard Model Higgs boson, provide a 
well motivated extension of the Standard Model (SM). In this work we show that, 
once new physics at heavier scales is added to stabilize the Higgs potential, such an 
extension of the SM can strengthen the first order electroweak phase transition and 
make the electroweak baryogenesis mechanism feasible. We propose a SM extension 
with TeV Higgsinos, Winos and Binos that satisfy the following properties: a) The 
electroweak phase transition is strong enough to avoid sphaleron erasure in the 
broken phase for values of the Higgs mass mn ~ 300 GeV; b) It provides large 
CP-violating currents that lead to the observed baryon asymmetry of the Universe 
for natural values of the CP-violating phase; c) It also provides a natural Dark 
Matter candidate that can reproduce the observed dark matter density; d) It is 
consistent with electroweak precision measurements; e) It may arise from a softly 
broken supersymmetric theory with an extra (asymptotically free) gauge sector; f) 
It may be tested by electron electric dipole moment experiments in the near future. 



1 Introduction 



In spite of all our recent progress in the understanding of physics at the electroweak scale, 
the source of Dark Matter (DM) and the origin of the matter-antimatter asymmetry pQ 
(BAU) still remain unclear. It is today well understood that the solution to either of these 
problems requires physics beyond the Standard Model (SM). 

On the one hand if the Standard Model baryon number violating interactions [3 El El 
cease to be in equilibrium in the bubbles of the broken phase, the matter- antimatter 
asymmetry may be generated at the electroweak phase transition via the mechanism of 
electroweak baryogenesis E| . For such a mechanism to be realized in nature a strongly 
first order electroweak phase transition is required. However the phase transition in 
the Standard Model for values of the Higgs mass consistent with the LEP bounds is 
a crossover [Zj and hence any baryon asymmetry generated at the weak scale would be 
erased. Moreover, the sources of CP- violation are insufficient to generate a baryon number 
consistent with the one observed in nature [8J. 

Electroweak Baryogenesis remains nevertheless as an attractive possibility in models 
of physics beyond the Standard Model at the weak scale. It has been shown that the mini- 
mal supersymmetric standard model (MSSM) is consistent with this mechanism provided 
one of the superpartners of the top quark is lighter than the top quark and the Higgs 
boson is lighter than ~120 GeV [H] JE]- This mechanism also demands the presence of 
charginos and neutralinos at the weak scale, which provide the necessary CP-violating 
sources [H| [HI] and also a natural Dark Matter candidate. Electroweak Baryogenesis 
may also be realized in the next-to minimal supersymmetric extension of the Standard 
Model (NMSSM) where some of the MSSM constraints can be relaxed. In particular there 
are modifications to the tree-level effective potential [32] that may ensure a strongly first 
order phase transition [33] [HE] without a light stop. 

On the other hand the Standard Model does not provide any natural source for the 
observed Dark Matter. Neutrinos are too light to give any sizeable contribution and 
there is no evidence of the possible existence of sufficient Jupiter-like, baryonic objects. 
Moreover, recent WMAP measurements |SH|-j40!] exclude the presence of a significant 
baryonic contribution to the observed Dark Matter density. The natural candidates for 
the source of Dark Matter are new, neutral, stable, weakly interacting particles with 
masses of the order of the weak scale. These particles lead naturally to a relic density 
of the order of the critical density and appear in many models beyond the Standard 
Model. In particular they are present in models of softly broken supersymmetry at the 
TeV scale. The lightest supersymmetric particle in these models tends to be neutral and 
its stability is ensured by a Parity symmetry, Rp, which also ensures the proton stability. 
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Acceptable values of the Dark Matter density and a successful realization of the mechanism 
of electroweak baryogenesis may be simultaneously obtained in minimal supersymmetric 
models, in certain phenomenologically interesting regions of the parameter space jH], [SB]- 
In the previous models, as well as in all successfully considered scenarios, the strength- 
ening of the phase transition proceeds from the existence of new, extra scalars in the 
theory, while the CP-violating and Dark Matter sources proceed from new fermion fields. 
The common lore from all previous works was that the presence of extra bosons was a 
necessary requirement to induce a strong enough first order electroweak phase transition. 
In fact it is currently understood that bosons coupled to the Higgs field <fr with coupling 
h favor a first order phase transition: they create a cubic term in the Higgs effective 
potential ~ (h<f)) 3 either at the tree level, as in the NMSSM, or by its contribution to 
the one-loop thermal effective potential ~ (h<f)) 3 T. On the other hand fermions do not 
give rise to any cubic term in the high temperature expansion in powers of h<p/T of the 
thermal integrals and hence they were neither believed to give rise to a barrier between 
the symmetric = and broken 7^ phases nor to trigger a first order phase transi- 
tion. In this paper we will prove that while the latter statement remains true for weakly 
coupled fermions h <C 1 it is not for strongly coupled (but still perturbative) ones h £ 1 
that can indeed induce a strongly first order phase transition consistent with electroweak 
baryogenesis. 

In this work we shall first show that in Standard Model extensions with extra fermions 
strongly coupled to the Higgs field the first order phase transition may be sufficiently 
strengthened in order to avoid erasure of the baryon asymmetry in the broken phase. We 
shall analyze in detail a simple model, which can be considered as a particular realiza- 
tion of split supersymmetry [32], where the standard supersymmetric relations between 
the Yukawa and gauge couplings are not fulfilled. We shall stress, however, that in 
such a model, the physical vacuum becomes unstable and therefore the strength of the 
electroweak phase transition may not be properly defined without an ultraviolet (UV) 
completion of this model, that includes the presence of heavier, stabilizing fields. An 
example of such fields may be provided by softly broken supersymmetry, although other 
extensions are possible. 

We shall show that this low-energy effective theory, with Higgsinos and gauginos 
strongly coupled to the Higgs, may arise from a soft supersymmetry breaking model, 
based on a gauge extension of the Standard Model gauge group, with new (asymptoti- 
cally free) gauge interactions that become strong at the TeV scale, and are responsible 
for the strong Yukawa couplings of Higgsinos and gauginos to the SM Higgs field. This 
gauge extension of the MSSM provides a UV completion of the model analyzed in this 
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paper and allows for large Higgs masses. 

The article is organized as follows. In section El we present the general ideas leading 
to the strengthening of the phase transition in the presence of strongly coupled fermions, 
and the need for the presence of stabilizing fields. In section El the phase transition for 
the Standard Model extension containing Higgsinos, Winos and Binos strongly coupled 
to the Higgs field is worked out in detail. We show that a strong enough first order phase 
transition can be accommodated even with a heavy Higgs, mn ~ 300 GeV. In section E] 
the CP-violating currents induced by the charginos that lead to the observed baryon 
asymmetry of the universe for natural values of CP-violating phases are presented. We 
show that in order to reproduce the WMAP results on the BAU the CP-violating phases 
must be 0(1) if all squarks are heavy enough to be decoupled from the thermal bath. 
Otherwise if some light squark (e.g. the right-handed stop) remains in the spectrum, 
the CP- violating phases can be as small as a few times 10 -3 . In section El the two-loop 
contributions to the electron electric dipole moment from the charginos and neutralinos in 
our model (in the absence of light squarks) are evaluated, assuming that all relevant CP- 
violating effects are associated with the new fermions. For the values of the parameters 
satisfying all other requirements the generated electric dipole moment is below the present 
experimental bound, although the model may be tested in the future if experimental 
bounds improve by a few orders of magnitude. In section El we discuss the Dark Matter 
constraints in our scenario, whereas compatibility of the strongly coupled fermions with 
electroweak precision measurements is considered in section El A natural region in the 
space of parameters is found where all requirements are fulfilled. A discussion on a possible 
UV completion of the model is presented in section El Finally we reserve section El for our 
conclusions. 



The finite temperature effective potential of the Higgs field is, by definition, the free- 
energy associated with 0. The one-loop, finite temperature contribution to the free-energy 
density is given by 



where is the number of degrees of freedom (d.o.f.) of the particle, mj(0) is the Higgs- 
dependent particle mass, and 



2 Phase Transition and TeV Fermions 




(2.1) 




(2.2) 
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where I-(x) [I + (x)] stands for the contribution from bosons [fermions]. In this section 
we will consider for simplicity masses of the form m 2 (0) = fi 2 + h 2 <j) 2 , where h is the 
Yukawa coupling and fi is a constant mass parameter 1 . For large masses (m/T ^> 1) the 
functions I T are exponentially small, which means that heavy species are decoupled from 
the thermal plasma. For small masses (m/T <C 1) I T can be expanded in a power series 
of m/T. This is the case of the minimal standard model. In fact adding this expansion to 
the zero-temperature effective potential in the SM one obtains the well known expression 
for the free-energy density 

Jsm (<P, T) = -^9*T 4 + V SM (0, T) , (2.3) 

where the first term comes from the entropy density contribution of relativistic particles, 
with the number of effectively light species in the plasma [13] (g* — 107 for the SM), 
and the second term is the field dependent effective potential 

V SM (0, T) — D (T 2 - T 2 ) <f - ET(j) 3 + ^<p\ (2.4) 

In the SM the parameters of Eq. (J2.4)) are well known and given by jUj 

D = ^{2m 2 w + m 2 z + 2m 2 ), E = JL (2m 3 w + m|) , 
T 2 o = B = ( 2m 4, + m 4 _ 4m 4) j 

()7T"r 1 
lr, /7„ — 3 

2w 2 ' 



A T = A - - - , 2m^ In + m| In — £j - 4m? In 



— — - + m 7 In — — r — 4m. In — — - 

a B T 2 z a B T 2 1 a F T 2 J 



m 



A = ^f, lna B = 3.91, lna F = 1.14. (2.5) 



where the vacuum expectation value of the Higgs at zero temperature is normalized to 
v ~ 246 GeV. 

A free energy of the form of the SM one, Eq. (|2.3jl . leads to a first order phase transition 
at a critical temperature given by 

T c = - T ° (2.6) 
T^EVA^D 



The value of the Higgs field at the minimum of the potential is 

3ET 



4>m (T) 



2A 



T 



(2.7) 



Mass eigenvalues are in general more complicated functions of </> as it will be the case in the model 
considered in section 
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and the order parameter at the critical temperature is given by (fi c /T c = <f) m (T c ) /T c — 



Notice that the parameter E is very small because only the (transverse) gauge bosons 
contribute to it and, as a consequence, <p c is much smaller than T c (unless A ^ 2E) and 
the phase transition is very weakly first-order 2 . Moreover, for physical values of the 
Higgs mass, the small value of <fi c /T c causes perturbation theory to break down and only 
non-perturbative calculations become reliable. To overcome this problem SM extensions 
containing extra bosons strongly coupled to the Higgs sector have been considered in the 
literature 0. In general these bosons would contribute to the parameter E and would 
strengthen the first order phase transition. In this paper we will prove that a similar 
effect is produced by fermions strongly coupled to the Higgs sector, even if in the high 
temperature regime they do not contribute to the cubic term in the effective potential. 

In general we will consider particle species that contribute to the one-loop effective 
potential as in Eq. (|2.1|) . not light enough for the validity of the power expansion of I T 
but not necessarily so heavy as to consider them to be decoupled in the typical range 
of temperatures of the electroweak phase transition. The ^-dependent part of the free 
energy density would then be given by 



F(cf>, T) = Jsm (0, T) ± ^(^M) + t4 J2 9ih K (4>) /T] /2tt 2 , (2.8) 



where the first term is given by Eq. ()2.3|) . the last one is the finite temperature contribution 
of the new, heavy particles, V(m 2 ) is the zero temperature contribution, and the plus and 
minus signs in front of V(m 2 ) correspond to bosons and fermions, respectively. The 
zero-temperature one- loop effective potential V(m 2 (4>)) is given by 



where P{4>) is a polynomial in that contains quadratic and quartic terms with coefficients 
that depend on the renormalization conditions [13]. By imposing the renormalization 
conditions already used in the SM, Eq. ()2.4j) . in particular that the tree-level values of 
the minimum and Higgs mass are not shifted by radiative corrections, i.e. 



2E/\ Tc . 



V{m 2 {(t))) 



1 



m 4 (^)logm 2 (0) + P{4>) 



(2.9) 



64vr 2 



dV 



0. 



d 2 V 





dcj) 2 



we obtain 




(2.11) 



2 For E = the transition becomes second order in the one-loop approximation. 
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with 



a 





UJijJ 


-3 
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+ uj' 2 + ujuj" ) logcu - - — + -uj' 2 + -ujuj" 



64vr 2 \ \ v ) & 2 v 2 2 

128ttV 



/5 = S 2 - u;' 2 - ujuj'') \oguj + — -3uj' 2 - ujuj''} (2.12) 



, and so on. For the case 

<f>=v 



where we are using the notation: uj = m 2 (v), uj' = dm d ^ 
where the <fi dependence of the mass eigenvalue is m 2 (0) = fi 2 + h 2 <ft 2 the potential (12. 9|) 
has the familiar expression 



V{m 2 (<P)) 1 



64^ 



m\<P) ( log ( 2m 2 (<p)m 2 (v) 



m 2 (v) J 2 



(2.13) 



Let us first stress that, unless the Higgs is heavy, strongly coupled fermions may create 
a problem of vacuum stability at scales close to the electroweak scale. This can be easily 
understood from the fact that the tree-level quartic coupling, defined as the coefficient of 
the quartic term in the effective potential, is given by m 2 H /8v 2 and the radiative corrections 
are proportional to the fourth power of the Yukawa coupling, h, and to the number of 
degrees of freedom. As we shall demonstrate, a relevant effect may only be obtained for 
a value of the number of degrees of freedom times /i 4 larger than O(10). For such values 
vacuum stability occurs at scales of the order of TeV and therefore the presence of new, 
stabilizing fields is necessary in order to define a consistent low-energy effective theory. 

An efficient way of stabilizing the potential in the presence of strongly coupled fermions 
is to assume the presence of heavy bosonic degrees of freedom with similar couplings and 
number of degrees of freedom. For simplicity, let us here assume that the fermions have a 
dispersion relation m 2 (<f)) = fi 2 + h 2 <f) 2 , and a number of degrees of freedom g, and there 
are bosonic, stabilizing fields with a dispersion relation m 2 s ((p) = /i| + h 2 (p 2 with equal 
number of degrees of freedom. Then, taking into account only the radiative corrections 
associated with these heavy fields, the maximum value of fi$ consistent with vacuum 
stability may be obtained from the condition of a positive quartic coupling at scales much 
larger than v , and it is given by 



^" 6XP {jh^J mf(v) ' k V (2 ' 14) 

Observe that for heavy Higgs bosons and/or weakly coupled fermions fis becomes much 
larger than the weak scale. If, however, h takes large values and m# becomes light then 
fis approaches fif and the effect of stabilizing fields is to cancel the zero temperature 
contribution of fermions plus giving additional finite temperature contributions. 
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In order to get an understanding of the effects to be expected by the presence of these 
new particles, let us first consider a fermion particle with a mass, m(4>) = h<p, much larger 
than T c for <fi ~ <p c . If this effect is obtained for large values of the Yukawa coupling, then 
the sum of the effects of the fermions and the stabilizing fields at zero temperature is 
small and, in this extreme case that maximizes the contribution of fermions (and that 
of the stabilizing fields) to the phase transition, we can ignore the zero temperature 
contributions. Then, in the symmetric phase the species is light and contributes to g*, 
but in the broken-symmetry phase it is heavy and approximately decouples from the 
thermal plasma. Usually such a decoupling species would transfer its entropy to the 
thermal bath, causing a temperature rise. During a phase transition we would naively 
expect that the effect of such a reheating is to delay the appearance of the true vacuum, 
to decrease the critical temperature and subsequently to increase the value of (f) c /T c . 

More quantitatively, at constant g* the critical temperature is given by the condition 
V {4>m (T' c ) , T c ') = 0, corresponding to degenerate minima of V 3 . As we will show below, 
this condition changes if the number of light degrees of freedom is different in the two 
phases, Ag* = g^ S y mme t r ic — 9* broken ^ ^- ^ n our exam pl e > t ne decoupling particle 
has a mass m(4>) = h<f), and its contribution to the free-energy density in the broken 
phase vanishes, while in the symmetric phase is equal to — ^-A^T 4 , where Ag* is its 
number of degrees of freedom. The condition of degenerate minima of J 7 , Eq. p. 8)1 . then 
gives [13 Ej 

^sm (<f> m (T c ) ,T C ) = -^A<7*T C 4 . (2.15) 

This condition is attained at a lower critical temperature, T c ' > T c . Moreover we have that 
<pc = 4>m (T c ) > <p' c and the phase transition is more strongly first-order. It was already 
noticed in Ref. jlH] that Ag* ^ in the context of the electroweak phase transition could 
be important for baryogenesis. 

We can now estimate A (<p c /T c ) = 4> C /T C — (p' c /T^ by noticing that the value of the 
effective potential (|2.4)1 at the minimum m can be written as 



V SM ((p m (T),T) A f<f>m(T)\ 3 f<f>' c <Pm(T) 



T 4 4 V T ) \TL T 



c 



(2.16) 



where we have made the approximation At — A. Eq. ()2.15jl implies that 



(2.17) 

T C J \Tj 45A 1 ' 

For the order parameter to increase from <fi' c /T' c = 2E/X <C 1 to the value 4> C /T C > 1, neces- 
sary to preserve the baryon asymmetry [5], we need Ag* > 45A/27T 2 ~ 0.25 (m^/115 GeV) 2 . 
3 We are normalizing the total effective potential as 1^(0, T) = 
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So it seems that particles with few d.o.f. will produce an effect. Conversely given Ag*, 
the bound on the Higgs mass is relaxed from (rnn/v) 2 < AE to 

(m H /v) 2 < AE + 4vr 2 A^/45. (2.18) 

Notice that if the new value <p c is ~ T c , then the perturbative (one-loop) approach is fully 
justified even for small values of E. 

It should be noticed that Eqs. (|2.15jl - (|2.17jl give only an estimate of T c and 4> c . In fact, 
the effective number of light d.o.f. varies continuously from to — Ag*, as <fi goes from 
to <p c . Therefore we have a function g m (0) that contributes to T (4>, T), and the correct 
value of the vacuum expectation value (VEV) is obtained by minimizing the complete 
free energy, Eq. (|2.8J) . Once the dependence of g* on <f> is taken into account, we cannot 
use analytical approximations and we need to resort to numerical calculations. Observe 
that in the limit studied above fermions and bosons gave equally important contributions 
to the phase transition strength and the number of degrees of freedom is the sum of the 
one associated to the fermions and that of stabilizing fields. In the rest of this work we 
will consider bosons that are heavier than the fermions, and therefore lead to a smaller 
finite temperature contribution than in the example above. 

We will now consider adding a fermion particle with mass m 2 (0) = fi 2 + h 2 <ft 2 to the 
SM with a Higgs mass m# = 120 GeV, but we will retain the effects of the heavy particles 
in the broken phase. We consider only fermion species and their (heavier) stabilizing fields 
since, as explained above, the effect of bosons on the phase transition has been extensively 
studied in the literature 0, JZJ-jBI]- ^ n Fig. ^ we plot the number of degrees of freedom 
g that give <p c /T c = 1 as a function of the mass parameter /i, for different values of the 
Yukawa coupling h. The invariant mass of the stabilizing fields has been set to their 
maximum value consistent with vacuum stability, Eq. (J2.14j) . As anticipated, for fi ~ 0, 
and for large values of the Yukawa couplings, only a small number of degrees of freedom 
are necessary in order to obtain a strongly first order phase transition. 

A minimal Standard Model extension (i.e. the introduction of a single species) is 
possible with bosons but not with fermions, since the SM Higgs is an SU(2) doublet. In 
order to construct an invariant Yukawa Lagrangian the simplest possibilities are, either 
doublet and singlet fermions (as e.g. a generation of mirror leptons and/or quarks) or 
doublet and triplet fermions (as e.g. light Higgsinos and gauginos) remnant from (split) 
supersymmetry 42,. In the next section we will consider the latter possibility. 
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Figure 1: Curves of constant 4> C /T C = 1 and m# = 120 Gel 7 , /or a fermion with mass 
m 2 = yU 2 + /i 2 2 and g degrees of freedom. From top to bottom the curves correspond to 
h = 1.5, 2, 2.5, and 3. 



3 HlGGSINOS AND GAUGINOS: THE PHASE TRANSITION 



In this section we will consider the particular case of Higgsinos (Hip), Winos and Binos 
(W a , B) coupled to the SM Higgs doublet H with the Lagrangian 

C = (ji 2 cr a W a + h' 2 B^j H 2 + H T e {^—hi a a W a + h[ B^j Hi 

M 2 ^ra^ra , M l 6 6 , .. tjT , 



+ ~y W a W a + -yBB + flH^ effi + /i.e. 

where e = ia 2 and the Yukawa couplings hi t2 and ^ 2 are arbitrary 4 . 
The chargino mass matrix is 

' M 2 hi<j 
h 2 (j) /i 

and the squared mass matrix has eigenvalues 

A,. = = ( J Ml + hlcj) 2 ± A /M 2 + hi<j) 2 * 



where M± = |(M 2 ± /x), /i± = ± (/i x ± /i 2 ). The mass matrix for neutralinos is 



/ M 2 


-h 2 (j)/V2 
V MA/2 





M x 
/i 2 ^/v^ 
-h' l( f)/V2 



-MA/2 MA/2 \ 



-// / 



(3.1) 



(3.2) 



(3.3) 



(3.4) 



4 The matching at high scale with the MSSM couplings would be h% = gsin/3/y/2, hi = 5 cos (3/ '\/2, 
/i' 2 = g' sm/3/y/2, h[ — g' cos (3/V2, where 5 and are the SU(2) and Z7 (1) gauge couplings respectively, 
and the Higgs doublet is related to the MSSM Higgses by H = sin {3H 2 - cos^eff* The matching 
with the couplings of a possible UV completion of the model will be done in section [SJ For the moment 
we will just assume that there is such a UV completion and that it provides the necessary stabilizing 
fields as it was discussed in section 2. 
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The eigenvalues of this matrix are cumbersome, so we consider the particular case M\ = 
M2 = M, hi = hi = h, and h[ = h' 2 = hi . The eigenvalues of the squared mass matrix 
are thus 

Am = Li 2 , Xn2 = M\ (3.5) 

and 

A n± = (m+ ± sJm 2 _ + h 2 <f) 2 + 2h' 2 <j)^J , (3.6) 

In this case the chargino eigenvalues become very similar to those in Eq. (J3.6|) : in particular 
we have h_ = 0, h + = h, and M 2 — M in Eq. (|3.3|) . We can further simplify the problem 
by also setting h! = 0, since in this case Eqs. ()3.3|) and ()3.6|) become 

A± = (m+ ± + h 2 ^ . (3.7) 

that corresponds to 6 degrees of freedom (a Dirac spinor and a Majorana spinor) with 
squared mass A + , 6 with squared mass A_, and 2 (light) Majorana particles with masses \x 
and M. A total of 16 fermionic degrees of freedom out of which only 12 are coupled to the 
SM Higgs. Clearly, for hi — the Majorana particle with mass M is just a pure Bino state 
of mass Mi = M, which decouples from the other low-energy states and therefore plays 
no role in determining the strength of the electroweak phase transition. In the following, 
we shall concentrate on this particularly simple case. 

From Eq. ()3.7|) it is clear that all the results in this case will be symmetric under 
fi <-> M. The simplest limiting case is when M + = and M_ = M = —fi. In this case the 
eigenvalues are degenerate, A± = M 2 + h 2 <p 2 , with 12 degrees of freedom corresponding 
to 3 Dirac spinors, and the situation is identical to the simple example illustrated in 
the previous section. One expects that other limits will be less favorable for the phase 
transition. For instance if we take M_ = (i.e. M + = M = /i), the eigenvalues are 
{M±h(j>Y. This means that, unless hep > 2M, in the broken-symmetry phase half of the 
particles become heavier than in the symmetric phase but the other half become lighter. 
Therefore, unless M + is also small, less degrees of freedom than in the case M + = will 
contribute to this effect. 

In Fig.Elwe plot the free energy at different temperatures for a Yukawa coupling h = 2, 
a Higgs mass of 120 GeV and \i = —M ~ 200 GeV. We can explicitly see from the figure 
shape that there is a first order phase transition with an order parameter <f) c /T c ~ 1.75. 

In Fig. El we plot the ratio of the Higgs VEV to the temperature, evaluated at the 
critical temperature, as a function of the mass M for the case \x = —M and mn = 
120 GeV. As expected, the strength of the phase transition decreases with M. This 
illustrates the fact that for large M the particle is decoupled already in the symmetric 



10 
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4>/T 

Figure 2: The free energy at different temperatures for h = 2, nin = 120 GeV and 
fi = -M ~ 200 GeV. 

phase, hence the VEV (f) c has a smaller value, which corresponds in the M — > oo limit, to 
that of the electroweak phase transition in the Standard Model. 



1 . 75 
1 . 5 
1 .25 

C /T C 1 
. 75 

0.5 

.25 



100 200 300 400 

M [GeV] 

Figure 3: C /T C for uih = 120 GeV and h = 2 as a function of M = —\l in GeV. 

In Fig. 0] we plot the values of the Yukawa coupling h necessary to induce a strongly 
first order phase transition for the case M = — \x and tuh = 120 GeV. It is clear from the 
plot that for such a small value of the Higgs mass, a strengthening of the phase transition 
may only be achieved for h > 1.5. Let us stress that this lower bound may be weakened 
by assuming slighly smaller values of us- For instance, taking /z| to be 0.9 times the 
maximum value allowed in Eq. ()2.14j) is enough to ensure that values of h = 1.5 and 
masses \i of order of 100 GeV are allowed. 

We are interested in a strong enough phase transition for baryogenesis, which means 
that the Higgs VEV must be large enough for sphaleron processes to be suppressed in 
the broken-symmetry phase (in order to avoid the washout of the BAU after the phase 
transition). In the case of Fig. [21 this happens up to a value M = M max determined by 
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Figure 4: Yukawa coupling h necessary to get <p c /T c = 1 for ran = 120 GeV as a function 
ofM = —fjL in GeV. 

the condition <p c /T c ~ 1. Observe that the top-quark, with 12 degrees of freedom (similar 
to the chargino-neutralino system) and a Yukawa coupling ht — 0.7 in our normalization, 
is not able to generate such a strong first-order phase transition in the SM. 

It is well known that the value of the Higgs mass plays a prominent role in the strength 
of the phase transition in the Standard Model and extensions thereof. Up to now we have 
fixed it to a "minimal" value, m# = 120 GeV. However our mechanism of strengthening 
the phase transition by using strongly coupled fermions, although certainly sensitive to 
the value of the Higgs mass, permits to go to higher values. In Fig. El we plot, for fixed 
values of the Yukawa coupling, the values of M that give <p c /T c = 1 as a function of the 
Higgs mass. We can see that, as expected, there is an upper bound on the Higgs mass 
that depends on the Yukawa coupling (h). Moreover the upper value of the 

Higgs mass has an approximate linear behaviour as a function of the Yukawa coupling h 
as can be readily deduced from Fig. El Imposing perturbativity of the theory at the low 
scale sets a generic upper limit on the Yukawa coupling h < y/Air ~ 3.5 which from Fig. El 
yields a corresponding upper limit on the Higgs mass of m# < 400 GeV. In general, for 
large values of h, the requirement of perturbative consistency of the theory up to high 
energies may only be fulfilled by embedding this model into a more complete theory where 
couplings become asymptotically free (see section [SJ). For the particularly interesting UV 
completion proposed in section |H] the upper limit on the Yukawa coupling is, as we will 
see, more restricted, h < 2, which translates into the stronger Higgs mass upper bound 
m H < 175 GeV for fj, = -M of order 100 GeV. 

The model we have proposed contains Higgsinos Hi and electroweak gauginos W a , 
B, coupled to the SM Higgs through the Lagrangian (|3.1|) with Yukawa couplings hi, h\ 




12 



150 200 250 300 350 
m H [GeV] 



Figure 5: Contours of <p c /T c = 1 in the (M, run) -plane for h = 1.6, 2, 2.5, 3 and M = —fi. 
The vertical line corresponds to the experimental lower bound, for a SM-like Higgs, of 
m H = 115 GeV. 

and Majorana masses Mj. We have concentrated on the case h! i = 0, in which the Bino 
state decouples from the other low energy states. Similar results for the phase transition 
would be obtained in any model in which the Bino would decouple from the low energy 
theory, for instance if h[ <C hi and M\ is much larger than the weak scale. Although the 
Bino is absent at low energies, there is still a state with mass approximately equal to 
that is a candidate for Dark Matter. This will be analyzed in section El where we will 
show that a modest splitting between the Yukawa couplings hi and h 2 will be necessary to 
accomodate the observed DM energy density, although a large separation is not permitted 
by electroweak precision measurements. In this way introducing a small splitting between 
h 2 and h\ is not expected to modify substantially the previous results in this section. In 
fact in the approximation of neglecting terms of order h 2 _ in the diagonalization of the 
mass matrices one can prove that the mass eigenvalues still correspond to one Majorana 
spinor with mass /i, two Dirac spinors with squared masses given in Eq. (J3.3|) and two 
Majorana spinors with squared masses, 



which of course coincides with the degenerate result of Eq. ()3.7)1 in the limit of h\ = h 2 - 
The result of Eq. (J3.8)) also proves that the modification of the strength of the phase 
transition due to the non-degeneracy will be small at least in the validity region where 
Eq. (|3.8J) holds. More explicitly, in Fig. |H1 we show numerically the variation of the order 
parameter <p c /T c away of the degenerate point as a function of 5h = h_ for h + = 2 
and M = —fi = 50 GeV. We see that the phase transition is weakened by < 20% for 




(3.8) 



h_/h + < 0.3. 
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Figure 6: Values of 4> C /T C as a function of 5h, for a fixed value of the average value h + = 2, 
m H = 120 GeV and M = —fi = 50 GeV. 

Up to this point we have discussed the electroweak phase transition in the absence of 
CP-violating phases. However for the baryogenesis mechanism to work we need a non- 
vanishing CP-violating phase in the parameters of the theory that will trigger baryon 
number generation. As in the MSSM studies, we shall consider real Yukawa couplings 
(related to gauge couplings in the MSSM) and Majorana gaugino masses, and a complex 
mass parameter \x = \ii\e tLp . To conclude this section we would like to consider how the 
phase transition and in particular the order parameter C /T C vary with the phase (p. The 
result is presented in Fig.[7|where we plot C /T C as a function of <p for M 2 = M 1 = \fi\ = 50 




<p (rad) 

Figure 7: <j> c /T c for h = 2 and \fi\ — M — 50 GeV as a function of <p = argQu). 

GeV and Yukawa couplings h + = 2, h_ =0. The plot at (p — tt is consistent with the 
corresponding results that can be read off from Fig. El while the strength of the first order 
phase transition decreases by only a tiny amount for 0(1) CP violating phases, i.e. for 



14 



<p = 7T/2. 



4 CP-Violating Sources and Baryogenesis 

The chargino sector in the model presented in section El has a similar structure to the 
chargino sector in the minimal supersymmetric Standard Model. The only difference 
is that the couplings gsin P/y/2 and gcos(5/\/2 are replaced by arbitrary couplings /12 
and hi, respectively, as can be seen in the corresponding mass matrix ()3.2j) . As in the 
MSSM the CP-violating phase can have its origin, after field redefinitions, in the phase 
<p of the (complex) /i-parameter. A general method for computing the effects of CP- 
violating mass terms on particle distributions was introduced in Ref. leading to an 
efficient transport of CP-violating quantum numbers into the symmetric phase where 
weak sphalerons are active and can trigger electroweak baryogenesis for all bubble wall 
widths. The method was adapted to the MSSM by a number of papers [HI- Ell where a 
set of coupled differential equations, that include the effect of diffusion, particle number 
changing reactions and CP violating terms, were solved to find various particle number 
densities diffused from the bubble wall, where CP-violation takes place, to the symmetric 
phase where sphalerons are active. These methods can be adapted to the present model by 
just considering the particular structure of the chargino mass matrix given by Eq. (|3.2j) . 
We will further make the simplifying assumption that all CP violation resides in the 
fermionic sector. Otherwise there should be extra contributions to the CP violating 
currents from the bosonic (stabilizing) fields, although these contributions are expected to 
be suppressed with respect to the fermionic ones because the stabilizing fields are heavier 
than the fermions. So from this point of view our results will be mostly conservative. 

In particular we will follow the formalism of Refs. j^Zl EH] where a method was de- 
veloped to compute the CP-violating sources induced by the passage of a bubble wall 
in a system of fermions that interact in a way similar to the one described above, in an 
expansion of derivatives of the Higgs fields. The method allows for the computation of 
the currents in a resummation to all orders of the Higgs vacuum expectation value effects. 
It was found that there are two different CP-violating sources from the chargino sector 
which the total baryon asymmetry depends upon. The leading contribution is provided 
by 

e ij H i d fl H j =v 2 (T)d fl /3 (4.1) 

that is proportional to the variation of the angle (3 = arctan [v2(T)/vi(T)] at the wall of 
the expanding bubble. The source (|4.1|) has a resonant behaviour for M 2 = |/i| and it is 
the leading contribution in the MSSM. However, the Higgs sector of our model (which 
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contains just the SM Higgs doublet) can be considered as the uia —>■ oo limit of the MSSM 
Higgs sector (where is the pseudoscalar mass) in which case 8^(3 — *> and the source 
(14. lj) does not contribute to the diffusion equations. 

The second contribution to the baryon asymmetry is proportional to 

Hxd^Hz + i^tfx = v 2 cos(2/?)<y? + vd^v sin(2/3). (4.2) 

In the limit rriA — > oo only the second term in ()4.2j) survives. Moreover we will consider 
in this section hi — h% = h (i.e. tan (3 ~ 1) 5 in which case the only remaining source is 
proportional to d^v 2 . This region is, as it was proven in Ref. |3*T] . very insensitive to the 
resonance region relating M2 and \fi\ and it provides a very natural region of parameters 
where electroweak baryogenesis can hold. Although in the MSSM such a region provided 
a very tiny amount of baryon asymmetry, in the present model all effects are enhanced 
by the strong Yukawa couplings of the Higgs to charginos as we will see in this section. 
In the MSSM there is an additional suppression of the source, Eq. ()4.2|) . due to the large 
values of tan (3 necessary to fulfill the LEP bounds on the lightest CP-even Higgs boson 
mass jlH], tan/? of order 10. As stressed above in this model, instead, tan/3 ~ 1. 

Following the formalism of Refs. |2Z1 EI] the solution of the diffusion equations, in the 
limit where the strong sphaleron (r ss ) and Yukawa processes (Ty) are fast enough, provide 
quark number density for third generation doublets uq and singlets as functions of the 
number density for the Higgs doublet coupled to the top quark n#, as 

/ \ k Q {9k T - k B ) ( 1 1 

k H {k B + 9k Q + 9k T ) \T SS T y 

( ^ fc T (9/c Q + 2k B ) ( 1 1 

k H {k B + 9k Q + 9k T ) \ Y ss T Y , 



where ki are statistical factors [SU] 



3 (" p 



2 



k B (m 2 ) = — / dp " — (4.4) 

2n Jo sinh 2 (^{p 2 + m 2 )/T 2 /2 

k F (m 2 ) = —- / dp g r (4.5) 



2vr 2 



cosh 2 v/(p 2 + m 2 )/T 2 /2 



that satisfy the condition kp(0) = 1 (k B (0) = 2) for Weyl fermions (complex bosons). 

In turn the density uh{z) 6 is obtained from the diffusion equations as a function of the 

5 This choice will be motivated by the contribution of charginos and neutralinos to the electroweak 
parameter T as we will see in section 

6 The spatial coordinate z is transverse to the bubble wall and we are neglecting the bubble curvature. 
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particle number changing rates, CP- violating sources and diffusion constants, as explained 
in Ref. [3T], yielding from Eq. ()4.3|) the quark number densities tiq^t- 

In order to evaluate the baryon asymmetry generated in the broken phase n B we first 
need to compute the density of left-handed quarks and leptons, til, in front of the bubble 
wall in the symmetric phase. These chiral densities bias weak sphalerons to produce a 
net baryon number. Considering particle species that participate in fast particle number 
changing transitions, and neglecting all Yukawa couplings except those corresponding to 
the top quark, only quark doublets do contribute to n L . Then assuming that all quarks 
have nearly the same diffusion constant it turns out that uq 1 = nq 2 = 2(nq + tit) 
and therefore from Eq. (J4.3)) 



n L (z) = 5n Q (z) + An T (z) = An H {z) + O 



1 1 



5 k Q k B + 8 k T k B - 9 k Q k T 

k H (k B + 9k Q + 9k T ) ' 1 " ' 

It turns out that the baryon asymmetry can be written as [2] 

n F T ws f° , , x ( 5n F T ws z\ 
n B = / dz n L {z) exp (4.7) 

where n F = 3 is the number of families, v w the bubble wall velocity and T ws = Qua^T, 
where k ~ 20 jl], is the weak sphaleron rate. 

For the model presented in this paper, where squarks and the non-SM Higgs bosons 
are superheavy, ttiq, mr, rn B ^> T c ~ 100 GeV {k B = kr — 3, Uq ~ 6 and kn ~ 8), it 
turns out that the coefficient in Eq. (|4.6J) is A ~ 0. This SM suppression was already 
pointed out by Giudice and Shaposhnikov (30] and consequently the baryon asymmetry 
n B in our model is produced by sub- leading effects. Assuming T Y 3> T ss we can go beyond 
the approximation of Eq. f)4.3|) and work out corrections of 0(1/T SS ). This was done in 
Ref. [T7] leading to an 0(1/T SS ) correction to n L (z), A ss n L (z), in our model as 

1 lZ 1 ss 

where D q ~ 6/T is the quark diffusion constant and the strong sphaleron rate is given 
by T ss = 6/t'|a^T, where k' is an order one parameter [TTj . When ()4.8|) is inserted into 
f!4.7|) it produces the baryon asymmetry generated by the sub-leading 0(1/T SS ) effects. 
We have numerically checked that this correction is insufficient to generate the observed 
baryon asymmetry of the universe. 

Another, more important, correction that can lead to a non-zero value of the baryon 
asymmetry to leading order in T ss are Yukawa and gauge radiative corrections to statistical 
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coefficients ki (or equivalently thermal masses) [50J. Expanding Eq. ()4.5|) in a power series 
of m 2 /T 2 we can write 

3 m 2 

k F (m 2 ) = 1 - |_ + 0(m VT 4 ) . (4.9) 

Keeping only the strong gauge (g s ) and top Yukawa {h t ) couplings one obtains in our 
model the statistical coefficients, 

k T = 3(1 + A s + A y ) 
k B = 3(1 + A.) 

k Q = 6(l + A s + iA y ) (4.10) 



with 



9i a 3^ 



, > ^y = ~id ( 4 - n ) 

where ^ ~ l/v^ is the top quark Yukawa coupling and g s ~ 1.2 the strong gauge 
coupling. Numerically, | A s | ~ 7 x 10~ 2 is more important than |Ay| ~ 2 x 10~ 2 but since 
the strong correction to all quarks is universal it cancels in the contribution to the baryon 
asymmetry. In fact to linear order in A, one can write the density n^z) in Eq. ()4.fi)l as 

n L (z) = -^-A Y n H (z) (4.12) 
lo 

and the proportionality coefficient turns out to be A ~ 4 x 10~ 3 . This is the reduction 

factor we get from such a sub- leading effect. The numerical calculation of the baryon- 

to-entropy ratio rj is presented in Fig. |H] (lower solid line) where we plot its ratio to the 

experimentally determined value 7/bbn = (8.7 ± 0.3) x 10 -11 [SS] as a function of the 

Yukawa coupling h and we have fixed the CP-violating phase sin (p = 1. In fact the phase 

sin (p that would be required for fixing rj = ^bbn is given by the inverse value plotted in 

Fig. |S1 We have chosen the case /i = — M 2 exp(^), and small values of M 2 = |/i| ~ 50 

GeV, where the phase transition is favoured, and typical values of the bubble width 

and velocity 7 . Since the computation of the baryon asymmetry has been performed by 

ignoring corrections of order one, the main conclusion one can extract from the results of 

Fig. |H1 is that CP- violating phases such that sin (p is of order one are necessary to obtain 

a value of the baryon asymmetry consistent with the experimentally determined values, 

for any value of h > 1.5. 

7 A general feature of first-order phase transitions is that the release of latent heat causes a slow- 
down of bubble expansion |48j . The electroweak bubble- wall velocity thus decreases during the phase 
transition [3] from its initial value v w ~ 1CP 1 52] given by the friction of the plasma. Calculating 
the exact value of v w is out of the scope of this paper. However, as noticed in Ref. the effect of 
the velocity variation on the BAU is likely to be an 0(1) one effect and should not modify the main 
conclusions of this paper. 
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Figure 8: The ratio ^/t^bbn as o, function of the Yukawa coupling h for fi = —M 2 exp(iip), 
M 2 = 50 GeV, maximal CP-violating phase, simp = 1, and bubble parameters = 10/T C , 
Vu = 0.1. Left-handed squarks and right-handed sbottoms are heavy (in the few TeV 
range). The lower (upper) solid line corresponds to heavy (light) right-handed stops, £ 
1 TeV (tut — 100 GeV). Dashed line corresponds to right-handed stops with ~ 500 
GeV. 

The amount of generated baryon number density can be increased if some squark is 
light enough to be in equilibrium with the thermal bath during the phase transition, 
in which case the SM suppression is avoided. The typical case that was considered in 
previous studies is that of a light right-handed stop P that corresponds to values of 
the supersymmetry breaking soft masses tuq, itlb ^> T c and ttit ~ T c . In that case 
the statistical coefficients are given by kQ ~ 6, kx — 9, ks — 3 and ku — 8 and the 
coefficient A in Eq. (14. 6|) does not vanish to leading order in 0(1/T SS ). In fact it is 
given by A ~ 1/6 and this produces an enhancement factor of (9(100) with respect to 
the case where all squarks are superheavy. This enhancement factor produces larger 
values of 77 (and so smaller phases are allowed) as can be seen in Fig. |H] upper solid line. 
Now fixing rj = i]bbn can be consistent with phases simp ~ 10 -2 . Notice finally that a 
similar enhancement would also appear if other squark species (i.e. right-handed sbottom 
or left-handed doublet) are light; this effect is not particularly linked to the lightness of 
right-handed stops. 

As mentioned above, the upper solid line in Fig. |H] corresponds to the extreme case 
where there are no extra bosons in the low energy spectrum, or equivalently stop masses in 



19 



the TeV range or larger 8 . Of course there can be intermediate situations where the stop is 
heavy but still does not fully decouple from the thermal plasma. In this case it contributes 
to the statistical factor kx with some small value that can significantly contribute to the 
A-parameter and departure its value from zero. For instance if mx = 5T C ~ 500 GeV, 
its contribution to the statistical factor k^ as given from ()4.4|) is ~ 0.24 which produces 
in (14. 6|) a value A ~ 1.2 x 10~ 2 and enhances the value of r\ from its value with fully 
decoupled right-handed stops. The corresponding value of r\ is plotted in Fig. |H1 (dashed 
line). We can see that the CP- violating phase for 7] = t?bbn i s now sin<y9 ~ 0.1. A similar 
effect would be produced by other not-so-heavy third generation squarks. 

5 Electron Electric Dipole Moment 

In the absence of light squarks, baryon number generation demands the presence of large 
phases in the chargino and neutralino sectors. In the previous section, we assumed gaugino 
masses and Yukawa couplings to be real, and therefore the relevant phase is the one of 
the \i parameter. Although one-loop corrections to the electron electric dipole moments 
are suppressed in the limit of heavy fermions, as has been stressed in Ref . jSl] , two-loop 
contributions become relevant. In this section we will evaluate the two- loop contribution 
to the electron electric dipole moment from the fermion and Higgs sector 9 . 

In Fig. we plot the chargino contribution to the electron electric dipole moment [53] . 
for fi = —M 2 exp(i(p), with M 2 real and positive (M 2 = maximal CP- violating phase, 
siny? = 1, h_ = 0, and for h + = 1.5, m H = 120 GeV (solid line); h + = 2, m H = 120 GeV 
(dashed line); and h + = 2, tjih = 200 GeV (dot-dashed line). We have verified that the 
results vary only slightly for non-vanishing values of \h^/h + \ < 0.1, which, as we will show 
in sections El and are preferred by dark matter and precision electroweak constraints. 
As discussed above, for values of h + < 1.5, the electroweak phase transition is weakly first 
order and the generated baryon asymmetry is not preserved. For slightly larger values, 
h + < 1.6, small values of |/i| < 100 GeV and Higgs mass values smaller than 125 GeV 
are demanded in order to make the phase transition strongly first order. The present 
electric dipole moment bound, \d e \/(ecm) < 1.7 10~ 27 (HEj does not put a bound on this 
model for these values of h + and For larger values of h + and similar values of 
Fig. El shows that the electron electric dipole moment contributions become smaller and 

8 For instance for mx = 10 T c ~ 1 TeV the stop contribution to kx is ~ 3 x 10~ 3 and A ~ 1.8 x 10~ 4 , 
much smaller than the previously considered thermal effects. 

9 We will assume here that squarks and stabilizing bosons are heavy enough not to contribute ap- 
preciably to electric dipole moments. However when particular UV completions of this model will be 
considered this assumption should be re-checked and if the new contributions are relevant they should 
be added to the fermionic ones. 
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Figure 9: Electric dipole moment of the electron as a function of \fi\ for maximal CP- 
violating phase, siny? = 1, and /i_ = 0. The Higgs boson mass and Yukawa couplings are 
fixed to run = 120 GeV, h + = 1.5 (solid line); run = 120 GeV, h + = 2 (dashed line); and 
m H = 200 GeV, h + = 2 (dot-dashed line). 

in addition, as shown in Fig. |H1 smaller phases are demanded for the baryon number 
generation. Therefore, the electric dipole moment bounds become even weaker in this 
case. 

However the anticipated improvement of a few orders of magnitude in this quantity [27| 
will be sufficient to test this model even for values of h + = 2, provided the values of are 
in the range necessary to obtain a good dark matter source and avoid the LEP invisible 
constraints: For large values of Mi, (or for M 2 = Mi) and small values of the lightest 
neutralino mass is mainly a Higgsino with mass approximately equal to For h + > 1.5, 
the chargino and the second lightest neutralino are heavier than 200 GeV and therefore, 
the stronger experimental bound comes from the Z invisible width. Surprisingly, for 
h_ — 0, the tree- level coupling of the lightest neutralino to the Z vanishes and therefore 
the invisible width bounds become very weak. However as we shall analyze in section® 
assuming thermal production of the lightest neutralino, an acceptable relic density may 
only be obtained for a small but non-vanishing coupling to the Z. Quite generally, if 
the dark matter density is determined by the s-channel Z-annihilation cross section of a 
neutralino, an acceptable relic density may only be obtained for values of the neutralino 
mass larger than 35 GeV, what in this case implies > 35 GeV j^H]. For such a range of 
values of a bound on the electric dipole moment of about 10~ 29 e cm will be enough 
to test this model for h + < 2 and simp ~ 0(1), as required for baryogenesis in the absence 
of light squarks. 
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6 Dark Matter 



One of the most attractive features of the model presented above is that the particles 
that lead to a strengthening of the electroweak phase transition are the same as the 
ones leading to a generation of the baryon number at the weak scale. It would be most 
important if the same particles would also provide a good Dark Matter candidate. As 
stated in the introduction stable, neutral, weakly interacting particles lead naturally to a 
Dark Matter relic density of the order of the one present in nature. As we will see, under 
the assumption of an i?-Parity symmetry, the lightest neutralino of the model presented 
above becomes a good Dark Matter candidate. 

Let us work in the simplest case, in which the Bino mass Mi takes large values (Mi ^> 
M 2 ) and mixes only weakly with the Higgsino {h' 12 <C ^1,2)- In such a case, due to the 
strong Yukawa couplings hi and h 2 , the charginos and two of the neutralinos acquire 
masses of about h + v. The mass of the lightest neutralino is close to and the lightest 
neutralino is therefore an almost pure Higgsino state. 

Assuming that all squarks, stabilizing bosons and heavy Higgses, if present, are con- 
siderably heavier than the lightest Higgsino, the states which determine the neutralino 
annihilation cross section are the light SM-like Higgs boson and the Z-gauge boson. Due 
to the small coupling of the SM Higgs boson to quarks and leptons, the annihilation cross 
section via s-channel Higgs boson production is very small, unless the neutralino mass 
is very close to m h /2, a quite unnatural possibility that we shall discard for the aim of 
this work. Hence the annihilation cross section is governed by the coupling of the lightest 
neutralino to the Z-gauge boson. 

The coupling of a neutralino state to the Z-gauge boson is proportional to the differ- 
ence of the square of the components of the neutralino into the two weak Higgsino 
states Hi, 

^oc(|A^| 2 -|^ 2 | 2 ) (6.1) 

This difference vanishes for values of hi = h 2 , and increases for increasing values of h_. 
Considering small differences between the values of hi and h 2 , the lightest neutralino, 
with mass approximately equal to is given by 

X ~ jh H 2 + / 2 Hi (6.2) 

and hence 

h\-h\ 

9 * z x WThl ■ (6 - 3) 

The annihilation cross section is proportional to the square of g^z and inversely propor- 
tional to the square of the difference between the lightest neutralino mass and the resonant 



22 



mass value, Mz/2. Therefore, the smaller the coupling g^z, the closer the modulus of the 
parameter /i should be to the resonant mass value. Hence, in order to get a value of the 
relic density consistent with the one determined by WMAP, there must be a correlation 
between the departure of |/x| from Mz/2 and the difference of the Yukawa couplings of 
the two Higgsinos to the Wino and Higgs field. 

The numerical estimates of the values of for a given value of h + have been obtained 
by computing the relic density, which is inversely proportional to the thermal average 
annihilation cross section, 



VLh 2 



1.07 x 10 9 GeV" 



M. 



pi 




dx WM gW 



(6.4) 



where x = m%/T, is the mass of the lightest neutralino particle and T the temperature 
of the Universe [SH]. The value of the variable x at the freeze-out temperature, Xf — 



m 



x/Tf, is given by the solution to the Eq. 



x f 



In 



0.038 (jg/gl?) m^M Pl (av) (x f ) 



x 



1/2 



(6.5) 



where g = 2 is the number of degrees of freedom of the neutralino, g* . is the total number 
of relativistic degrees of freedom at temperature Tf , and Mpi is the Planck mass. The 
thermal average of the annihilation cross section may be computed by standard methods 
and, for a particle of mass m^, is given by 



(av) 



! ds Js - 4m| W Ki(y/s/T) / 1Q m\ T K 2 (m ji /T), (6.6) 

J Am\ V ' I 



where s is the usual Mandelstam parameter, K\ and K% are modified Bessel functions, 
and the quantity W is defined to be 



W 



n 



d 3 pf 



(2 7 r)V 4 )(p 1 +p 2 -^p / ) \M\\ 
f 



(6.7) 



where \A4\ 2 is the squared matrix element averaged over initial states, and summed over 
final states. 



In Fig. we plot the required value of \fi\ in order to obtain the central value of 
the relic density consistent with the recent experimental results Qh 2 ~ ^wmap^ 2 = 
0.11 ± 0.01 59 as a function of h-. The average values of the Yukawa couplings have 
been fixed to h + = 1.5 (dashed lines) and h + = 2 (solid lines). Lower and upper lines 
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Figure 10: Values of the fi parameter leading to a value of the neutralino relic density 
consistent with the central value of the WMAP observation, as a function of h_, for 
H = —M2exp(iip) and M2 = |//|, assuming that the average value of the Yukawa couplings 
h + is fixed to the value h + = 2 (solid lines) and h + = 1.5 (dashed lines). Lower and upper 
lines for each set of h + values correspond to if = and if = tt/2, respectively. 

for each set of h + values correspond to two different values of the CP-violating phase, 
if = and if = tt/2 respectively. Larger values of \h_\ lead to larger values of We 
show the results only for small values of \h_\, since as we will discuss in the next section, 
large values of lead to unacceptably large corrections to the precision electroweak 
observables. A non-vanishing phase has only a significant impact on the relic density 
determination for values of h + ~ 1.5 and > 80 GeV ~ 0.2), for which the phase 
transition is no longer strongly first order and/or unacceptably large corrections to the 
precision electroweak data are generated. Otherwise, the variation induced by the CP- 
violating phases is of the order of (or smaller than) the present experimental uncertainty 
in the relic density Indeed, a one sigma variation of the relic density results for values of 

~ 80 GeV (60 GeV) would be obtained by varying by about 3 GeV (1.5 GeV). 

It is important to compare the results shown in Fig. El with those in Fig. For 
h + < 1.6 it follows from Fig. that in order to preserve a strongly first order phase 
transition one needs | | < 100 GeV and a Higgs mass smaller than 125 GeV. As seen in 
Fig. El such small values of are also consistent with Dark Matter constraints, provided 
that ^ 0.3. As emphasized before, and as we will show in detail in section 7, the 
requirement of consistency with the precision electroweak data further constraints the 
allowed parameter space. 

For h + = 2, the constraints coming from the requirement of a sufficiently strong phase 
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transition are much weaker. Values of and of the Higgs mass as large as 200 GeV are 
consistent with this constraint. This is a much larger region of values of the parameter 
\n\ than the one consistent with Dark Matter constraints, for the small values of that 
are required by precision electroweak data (see section 7). Therefore, also in this case the 
allowed region of parameters may only be determined once the analysis of the constraints 
coming from the consistency with precision electroweak data are evaluated. 

Finally notice that, for a fixed value of h+ in Fig. the region above (below) the 
line Q = W map corresponds to Q > ^wmap (^ < ^wmap)- Therefore while the region 
above the corresponding curve is excluded since it would predict too much DM density, 
the region below it is not excluded provided there is another candidate for Dark Matter 
in the theory. 

7 Electroweak precision measurements 

Heavy fermions, with large couplings to the Higgs field, may induce large corrections to 
the electroweak precision measurement parameters. Extra contributions may come from 
the stabilizing fields, but in this section we shall assume that the UV completion of the 
model is such that they are small. A very trivial (ad hoc) way of achieving this is if the 
stabilizing fields are a set of scalar singlet fields strongly coupled to the Higgs boson. For 
instance we can consider the case of N scalar (complex) singlets giving a contribution to 
the scalar effective potential given by 

V s = h 2 \S\ 2 \H\ 2 + fi 2 s \S\ 2 + X S \S\\ (7.1) 

with Ag > 0. This set of singlet fields contributes to the one-loop effective po- 
tential of the Higgs field as gy = 2N bosonic degrees of freedom with a mass squared 
m 2 ((f)) = fi 2 s + h 2 (f) 2 that corresponds to the typical case of stabilizing fields introduced in 
section|21 10 . As it is obvious this system of stabilizing (singlet) fields would not contribute 
to the electroweak precision observables or to the CP-violating observables analysed in 
the previous sections. In the following, we will just concentrate on the contribution to the 
electroweak precision measurement parameters from the fermionic sector of the theory. 

It is well known that if the fermion masses proceed from the usual contraction of a 
left-handed fermion doublet and a right-handed fermion singlet with the Higgs doublet, 
and if the mass difference of the fermion components of the doublet field is small, the 

10 Strictly speaking, for the stabilization of the effective potential, it is not necessary that the bosonic 
and fermionic number of degrees of freedom, gj, and g/, are equal and/or that the corresponding Yukawa 
couplings hb — hf, as we have been assuming in section |3 It is easy to see that a necessary stabilization 
condition is provided by gbh\ > gj /ji. 
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contribution of heavy fermions to the ^-parameter is about I/671". On the other hand the 
contribution to the T-parameter would depend on the size of the mass difference between 
the up and down fermions. Cancellation of anomalies requires the presence of at least 
two such new heavy fermion doublets and therefore the contribution to the ^-parameter 
tends to be large, S > 0.1. 

In the case under analysis, however, the symmetry breaking masses proceed from the 
coupling of the Higgsino doublets to the Higgs field and an SU (2) ^-triplet of Winos. 
Contrary to the standard case of heavy fermions, the contribution to the S-parameter 
becomes small in this case. For a given value of the average Yukawa coupling h + , the 
contribution to the T parameter, instead, becomes sizeable for large values of h_, while the 
contribution to the parameter U is an order of magnitude smaller than the contribution 



We shall work in the limit in which M\ is large. Thus the mixing of Binos with the 
Winos and Higgsinos becomes small and therefore the Binos decouple from the precision 
measurement analysis. For large values of the Yukawa couplings, the lightest neutralino 
has a mass close to |/i| and a coupling to the Z-gauge boson given by Eq. ()6.3|) . As 
explained in section El for h- = there are two Dirac charginos degenerate in mass with 
two Majorana neutralinos and the T parameter contribution vanishes. The mass difference 
between the neutralinos and charginos grows linearly with h_, as does the coupling of the 
lightest neutralino to the Z-gauge boson, and the T parameter grows quadratically with 
as shown in Fig. for values of h + = 2 (solid line) and h + = 1.5 (dashed line). 

Moderate contributions to the T-parameter are not in conflict with electroweak pre- 
cision measurements. Even in the absence of any other physics at the weak scale the 
corrections to the T-parameter coming from the neutralino and chargino sector may be 
largely compensated by the negative contribution induced by the presence of a heavy 
Higgs, which contributes to the S and the T parameters in a way proportional to the 
logarithm of its mass, 



where rnh ref is a reference Higgs mass value. 

The model under analysis falls therefore under the class of models which give a small 
contribution to the U parameter and sizeable contributions to the T parameter. Although 
the new physics gives only negligible contribution to S, a sizeable contribution to the S 
parameter may also be induced by a heavy Higgs boson. A fit to the precision electroweak 



to T. 




(7.2) 
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Figure 11: Chargino and neutralino contributions to the T -parameter as a function of h_, 
assuming that the average values of the Yukawa couplings has been fixed to h + = 2 ( solid 
line) and h + = 1.5 (dashed line) and considering M = —fi determined as a function of 
h_ from Figure [771 

data in this class of models has been done by the LEP electroweak working group [HO] . 
For a reference Higgs mass value of 150 GeV, they find 

S = 0.04 ±0.10 

T = 0.12 ±0.10 (7.3) 

with an 85% correlation between the two parameters. Taking this into account we show 
in Fig. [T2]the 68% C.L. (solid lines, hatched ellipses) and 95% C.L. (dashed lines) region 
of allowed values of the new physics contribution to the S and T parameters for values 
of the Higgs boson mass equal to 115 GeV, 200 GeV and 300 GeV, respectively Observe 
that, since we are presenting a fit to the allowed new physics contribution to the S and T 
parameters, for each value of the Higgs mass the origin of coordinates represents the SM 
value. 

From the results of Figs. E3 [TH an d El we see that for a fixed value of the average 
Yukawa coupling, h + , the requirement of an acceptable relic density and a good fit to 
the precision electroweak data implies an interesting correlation between values of h_, 
the parameter and the Higgs mass. The value of the Higgs mass is also bounded 
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Figure 12: 68 % (solid lines) and 95 % (dashed lines) C.L. allowed values of the S andT 
parameter based on the fit to the precision electroweak data in Ref. JEHj , for three different 
values of the Higgs mass: 115 GeV, 200 GeV and 300 GeV. 

from above from the requirement of a successful generation of the baryon asymmetry. For 
instance, for h + = 2, Fig. shows an upper bound on the Higgs mass of about 200 GeV 
and therefore from Fig. EH it follows that a good fit to the electroweak data may only be 
obtained for values of the new physics contribution to the parameter T < 0.27 (if the new 
physics contributions to the parameter S are small, as in the case under analysis). The 
upper bound on T translates into an upper bound on h_. From Fig. ^2 w e get that < 
0.18, and from Fig. El this leads to an acceptable relic density only for \fi\ < 65 GeV. 

For h + < 1.6 the corrections to the precision electroweak parameters are slightly 
smaller than for h + = 2. However, as emphasized before, from Fig. El we see that consis- 
tency with a strong electroweak phase transition can only be obtained for a Higgs mass 
value close to its present lower bound, rrih < 125 GeV, and small values of \fi\ < 100 GeV. 
From Fig. EE this implies that the new physis contribution to the parameter T < 0.15, and 
therefore from Fig.lTTlone obtains that < 0.14. Interestingly enough, from Fig.^J we 
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find that the same region of parameters is consistent with the Dark Matter relic density 
provided |/i| < 70 GeV. Observe that the allowed values of approximately coincide 
with the ones obtained for h+ = 2. 

8 Large Yukawa couplings in low energy super- 
symmetric THEORIES 

In this work we have analyzed the properties of a model with light charginos and neu- 
tralinos which are strongly coupled to the SM Higgs. Such large values of the Yukawa 
couplings do not arise in minimal supersymmetric extensions of the Standard Model, since 
these couplings are determined by supersymmetric relations to the weak couplings. It is 
therefore an important question to understand under which conditions such a low energy 
effective theory may be obtained. 

One possibility is to assume that, although the quantum numbers of the light particles 
are those of Higgsinos and gauginos of supersymmetric theories, the theory is not related 
to any supersymmetric theory at high energies. A more interesting possibility would be to 
consider this theory as a particular realization of split supersymmetry [12] , but where the 
particular relation between the Yukawa couplings and the gauge couplings has been broken 
by supersymmetry breaking interactions. One of the problems with this alternative is the 
one of vacuum stability. In addition, due to the strong Yukawa couplings, perturbative 
consistency is lost at scales much lower than the GUT scale and therefore perturbative 
unification of the gauge couplings cannot be achieved. 

In this section we shall show that this low energy effective theory may also arise in 
low energy supersymmetric models with extra, strongly coupled gauge sectors, if the scale 
of supersymmetry breaking is larger than the scale of spontaneous symmetry breakdown 
of the extended gauge sector to the Standard Model one. A supersymmetric extension 
provides the necessary stabilizing fields in a natural way. On the other hand in this exten- 
sion the strong Yukawa couplings are proportional to the strong gauge couplings, which 
become asymptotically free at high energies. This ultraviolet completion of the theory 
allows the preservation of perturbative consistency up to high energies and therefore the 
possibility of perturbative unification of gauge couplings. 

In order to illustrate this property let us consider the model analyzed in Ref . [EI] • The 
model is based on the low energy gauge group 

SU(3) C ® SU(2) 1 ® SU{2) 2 ® U(l) Y . (8.1) 

First and second generation left-handed quark and leptons transform in the fundamental 
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representation of S77(2)i, i.e. (2, 1) of S77(2)i ® SU(2) 2 , while the third generation left- 
handed quark and leptons and the two MSSM Higgs bosons transform in the fundamental 
representation of SU(2) 2 , (1,2). The model also includes a bifundamental (2,2) Higgs 
field E, as well as a singlet field S. The scalar component of the bifundamental takes a 
vacuum expectation value (E) — u ■ I, breaking S77(2)i ® SU{2) 2 — > SU{2) L . The gauge 
bosons 

W » = 9.W - 91 W» 
V9 2 i + 92 

remain massless after this symmetry breakdown and interact with an effective gauge 
coupling gw = 9\92l \/ 9\ + 9 2 - This model can be made consistent with gauge coupling 
unification by embedding the group SU(3) C ® S77(2)i g) SU{2) 2 <g> U{1)y into the grand 
unified product group SU(5) ® SU(5) broken by bi-fundamental field (5,5) diagonal 

vev's jnUES- 

In the model of Ref. [HI] there are extra fields transforming under SU (2)i but not 
under SU{2) 2 . With this particle content the coupling g 2 of SU{2) 2 is asymptotically 
free, but g± of SU(2) 1 is not. We shall work under the assumption that the coupling 
g 2 becomes strong at the scale u of spontaneous symmetry breaking of the symmetry 
SU {2)i ® SU (2) 2 , and therefore the Winos of SU (2) 2 interact strongly with the Higgs and 
Higgsinos. For our analysis here, the relevant term in the superpotential is W = M^EE, 
and the relevant supersymmetry breaking terms in the gaugino-Higgsino sector are the 
masses of the gauginos, M.^ corresponding to the two groups SU{2)i. In the absence of 
supersymmetry breaking terms the superpartners of would become the low energy 
Winos, interacting with Higgsinos and Higgs bosons with the weak coupling gw- In 
order to get a strongly interacting Wino-Higgsino-Higgs sector at low energies we need to 
decouple the bifundamental Higgsinos E, as well as the weakly interacting Wino W\ from 
the weak scale theory. This may be achieved by choosing the supersymmetry breaking 
parameters Aii and the Higgsino mass parameter to verify 

Mi, M s > g 2 u, M 2 < g 2 u (8.3) 

Supersymmetry is therefore broken before the scale of breakdown of SU{2)i ® SU{2) 2 
group to SU{2) l- The Wino of SU (2)i is only weakly coupled with the bidoublet Higgsinos 
and its large supersymmetry breaking mass ensures its decoupling from the low energy 
theory. 

For the parameters given above the low energy Wino has a mass 

M 2 ~ M 2 - %- (8.4) 
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and has a component on the strongly coupled Wino of SU "(2) 2 of order cos^s, with 



sin6» E ~ g 2 u/Mx. 



(8.5) 



The effective Yukawa couplings between the low energy Winos and the two Higgsinos are 
therefore given by 



where we have assumed that the CP-odd Higgs mass is larger than the weak scale, and 
therefore a single, SM-like, CP-even Higgs boson remains in the low-energy theory. As- 
suming a 2 = 92 ~ 1 we § e t that, in this realization of the low energy theory, the 
Yukawa couplings hi < 2. 

On the other hand, a very large supersymmetric mass Ms would also demand to be 
compensated in a precise way by a similarly large supersymmetry breaking mass in order 
to get the proper SU{2) 1 ® SU(2) 2 breaking scale. Lower values of M s would reduce 
the strongly coupled gaugino component of the low-energy Wino state and would require 
a fine-tuning between the two terms in Eq. ()8.4j) in order to obtain a small value of the 
effective low-enegy Wino mass. Therefore in this model, a moderate amount of fine-tuning 
is required to get consistency with phenomenological constraints and a strong electroweak 
phase transition 11 . 

We can now ask if the strongly coupled gauge bosons may serve as the stabilizing 
bosonic fields defined in section |21 Let us first stress that the particular extension of the 
MSSM presented above leads to extra contributions to the precision electroweak observ- 
ables. Small values of these extra contributions may only be obtained for values of g 2 u 
larger than a few TeV [SUESl- Since g 2 u acts as the bosonic mass [is defined in Eq. (j2.14|) . 
these particles can only act as stabilizing fields if the Higgs is heavier than the range con- 
sistent with a strongly first order phase transition. It would be interesting to investigate 
possible regions of parameter space in which cancelations between different contributions 
take place and consistency with data may be achieved for lighter gauge boson masses. 
Otherwise, additional fields would be necessary to stabilize the Higgs potential. 

A very simple possibility, that has already been pointed out in section is that the 
stabilizing fields are gauge singlets, harmless from the point of view of electroweak pre- 
cision measurements. These singlets should couple strongly to the Higgs sector, in order 
to stabilize the fermionic part of the zero-temperature effective potential, but should not 
contribute strongly to the Higgs quartic coupling since in that case the phase transition 
11 We thank D. Morrissey for helpful discussions on this point. 



hi ~ g 2 cos 9y, cos f3 / \/2 
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would become weakly first order. A simple example of a model producing the required 
effect is given by a singlet superfield P coupled to the Higgs field doublets as well as to a 
set of TV singlet fields Si, with superpotential 

W = \iS 2 P + X 2 PH 1 H 2 + -M S S 2 - l -M P P 2 + /i#i# 2 (8.7) 

where all couplings and the masses Ms and Mp are positive. The absence of a coupling 
of the singlet fields Si to the Higgs doublets as well as the appearence of only terms 
proportional to S 2 in the superpotential may be understood as a result of the invariance 
of the theory under a global O(N) symmetry. 

The supersymmetric masses of the singlet fields, Ms and Mp, are assumed to be 
much larger than the weak scale, suppressing the mixing of light singlet fermions with the 
standard gauginos and Higgsinos. We shall also assume that there are supersymmetry 
breaking effects in the bosonic S'-sector that prevent the possibility of integrating out 
the superfields S from the weak-scale theory and allow the bosonic ^-sector to be the 
stabilizing fields with mass given by (|2.14j) . Such supersymmetry breaking effects should 
also ensure the preservation of the modifications to the low-energy Higgs quartic coupling 
induced by the presence of the superfields S in the theory. 

Furthermore we will assume that the superfield P can be integrated out supersym- 
metrically so that for scales below Mp it gives rise to an effective superpotential as 

W cS = (AxS 2 + X 2 H 1 H 2 y + hl s S 2 + fiH 1 H 2 (8.8) 

The above superpotential gives rise to a coupling in the tree-level potential between the 
S sector and the Higgs sector, /i 2 |S'| 2 |iJ| 2 , where 

^~^A 2 A lS in2/3, (8.9) 

to an if-quartic coupling as |AA|if | 4 , where 

AA ~ -^A 2 sin 2/3, (8.10) 
M P K 1 

and to a self-interacting quartic coupling As|£| 4 fields 

X S ~ gL$. (8.11) 

Since we have concentrated in this work on the results obtained in the model with 
strongly interacting gauginos and Higgsinos for small values of h-, the value of tan/3 ~ 1, 
and therefore sin 2(3 ~ 1. For Ms ~ Mp and A2A1 > 1, one obtains values h\ > 1, 
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necessary for the singlets Si to stabilize the one-loop Higgs potential against the fermion 
contribution. On the other hand, for fi <C Mp, as required to obtain a good description 
of the dark matter relic density and precision electroweak observables, the contribution 
to the quartic coupling is not too large. This has implications for the Higgs boson mass. 
In fact, for values of tan/3 ~ 1, the D-term contribution to the Higgs mass is small and 
the Higgs mass may be raised above (or around) the MSSM values, due to the tree-level 
contribution 

m 2 H ~ 2A\v 2 + (loop - effects) + (D - term) . (8.12) 

Values of the Higgs mass of order 115-200 GeV may be naturally obtained for A|///Mp ~ 
0.1-0.3, that appear in this model for values of Mp of the order of the TeV scale and 
values of A2 somewhat larger than one. By choosing the appropriate value of N and the 
values of the couplings \, one can reach a situation similar to the one described in the 
previous sections. A detailed analysis of the parameter space in this particular model is 
outside the scope of the present paper. 

9 Conclusions 

In this article we have shown that heavy fermions with strong couplings to the Higgs fields 
may induce a strengthening of the electroweak phase transition and can also provide the 
proper CP-violating sources for the generation of baryogenesis. These heavy fermions, 
however, also induce for light Higgs bosons an instability of the Higgs potential at zero 
temperature and therefore require an ultraviolet completion of the theory to recover the 
consistency of the low-energy theory. In this work, we have assumed that the heavier, 
stabilizing fields have oposite statistics but similar couplings and number of degrees of 
freedom as the fermion fields. The above properties are then associated with the low 
energy theory consisting both of the heavy fermions and the heavier, stabilizing fields. 

We have illustrated this possibility by considering a model with TeV scale Higgsinos 
and gauginos that may lead to a sufficiently strong first order electroweak phase transition 
for values of the Higgs mass as large as 300 GeV. This is quite different from the results 
of the MSSM, in which a light stop is necessary, and the Higgs mass should be lower than 
~ 120 GeV to enhance the strength of the electroweak phase transition. Also at variance 
with the case of the MSSM is the fact that in this scenario the particles that induce a 
strong first order phase transition are the same ones responsible for the generation of the 
baryon asymmetry at the weak scale. 

This model preserves most of the properties of low energy supersymmetry, including 
a good Dark Matter candidate. Beyond the problem of vacuum stability, however, the 
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low-energy Yukawa couplings deviate from the ones obtained in minimal supersymmet- 
ric extensions of the SM, and their strength spoils the perturbative consistency of the 
theory at scales of about a few TeV. We have shown, however, that the model may be 
considered as the low energy effective description of a gauge extended supersymmetric 
standard model, in which the strength of the Yukawa couplings is related to the gauge 
couplings of an extended asymptotically free gauge sector, that becomes strongly inter- 
acting at TeV scales. This ultraviolet completion of the model solves the strong coupling 
problem and introduces new heavy particles that, due to their supersymmetric relations 
to the Higgsinos and gauginos, tend to stabilize the Higgs potential in a natural way. The 
stability of the potential, however, may only be achieved for values of the heavy particle 
masses smaller than about 1 TeV. These relatively small values of the heavy gauge boson 
masses lead to large corrections to the precision electroweak observables, unless cancella- 
tions between different contributions occur. Alternatively, we also have shown that the 
model can be completed with singlets that may stabilize the zero-temperature effective 
potential while providing Higgs masses consistent with present experimental bounds and 
with the required strongly first order phase transition. Finally, we have shown that this 
model may be tested by electron electric dipole moment experiments in the near future. 
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